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We study the order of the phase transition in the 3d U(l)+Higgs theory, which is the Ginzburg-Landau theory 
of superconductivity. We confirm that for small scalar self-coupling the transition is of first order. For large scalar 
self-coupling the transition ceases to be of first order, and a non- vanishing scalar mass suggests that the transition 
may even be of higher than second order. 



1. Introduction 

Scalar electrodynamics (i.e., the U(l)+Higgs 
theory) in three dimensions describes several dif- 
ferent physical systems. First, it is an effective 
high-temperature theory for 4d scalar electrody- 
namics (with or without fermions) Q. In par- 
ticular, the 3d theory can be used to study the 
finite temperature phase transition in the 4d the- 
ory for small 4d coupling constants. Second, it 
is directly the Ginzburg-Landau theory of super- 
conductivity. In this well, the properties 
of the phase transition between the normal and 
superconducting states are of much interest ® (tJ . 

In general, the properties of the phase tran- 
sition in the 3d U(l)+Higgs theory cannot be 
reliably studied in perturbation theory. This is 
due to the fact that the 3d coupling constants 
are dimensionful, resulting in infrared problems 
if there are nearly massless excitations (the 
true dimensionless expansion parameter is a cou- 
pling constant divided by a mass ~ e|/m_fr(T)). 
While the IR-problems are not quite as severe 
as in non-abelian theories containing gauge self- 
interactions, they nevertheless exist. In addition, 
non-perturbative topological defects play a role 
in U(l)+Higgs. 

Due to the difficulties mentioned, the 3d 
U(l)+Higgs theory should be studied non- 
perturbatively More specifically, some im- 

portant problems to be addressed are: (a) the ex- 
istence and order of a phase transition as a func- 
tion of the scalar self-coupling, (b) the character- 



istics of the phase transition: the critical temper- 
ature T c , latent heat, surface tension, correlation 
lengths in the two phases, and (c) the conver- 
gence of perturbation theory in the broken phase. 
We present here preliminary lattice results on (a); 
more complete results on (a) and (b) will be pre- 
sented elsewhere [||. (c) was addressed in ||. 

2. The action and the parameters 

The continuum U(l)+Higgs theory in 3d is de- 
fined by the action 



S 



<Px 



FijFij + iDi^iDi 



(1) 



where Di — di + ie^Ai. The scale of the theory is 
given by the dimensionful [GeV] gauge coupling 
e§, and the other two parameters are dimension- 
less, 



i(e§) 



A 3 

e 3 



The parameter x is roughly proportional to the 
ratio of the squares of the scalar and vector 
masses deep in the broken phase, x ~ m^/m^/, 
while y is related to temperature, y ~ (T— T c )/T c . 
Here T c is the tree level critical temperature. 

The Ginzburg-Landau theory of superconduc- 
tivity is defined just by eq. (|l]). Conventionally, 



A3 are denoted by a, b in that context, and x is 
represented by the Ginzburg parameter k = y/x. 
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If « < 1/V2 the superconductor is of type I, if 
k > l/v2 it is of type II. The phase transition 
is of first order for strongly type I superconduc- 
tors, and for strongly type II superconductors it 
is usually assumed to be of second order. 



• x=0.0463 

□ x=2 




3. Discretization and simulations 

The action in eq. (|l|) can be discretized in the 
usual way. We use the compact formulation for 
the gauge field, so that the 3d lattice action is 

S = G ^Re (1 - Up) 

p 

x,i 

+E«+Efe[ 1 -^(^wi 2 - ( 2 ) 

X X 

Here Ui(x) is the compact link variable, Up is 
the product of link variables around a plaquette 
P and <f> is a complex scalar field located on sites. 
It would also be possible to use a non-compact 
formulation Hf|. 

The three coupling constants in the lattice ac- 
tion are related to continuum parameters through 
a constant physics curve, which can be calcu- 
lated using lattice perturbation theory. Due to 
the fact that the couplings are dimensionful, the 
relation can be found exactly with a 2-loop cal- 
culation |ic[| . The relevant relations are 
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x (In 6/3 G + 0.09) + 25.5 + 4.6.x 



The continuum limit is hence at (Pg> Ph, Pr) = 
(oo,l/3,0). Note that increasing Ph means de- 
creasing y. 

The simulations were done using Cray C94 at 
the Center for Scientific Computing in Helsinki. 
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Figure 1. The maximum of the susceptibility \ as 
a function of volume. The straight line indicates 
the linear behaviour characteristic of first order 
transitions. The data point corresponding to the 
largest volume for x = 0.0463 was obtained by 
reweighting the data. 

4. The order of the transition 

For small x, the phase transition in the 3d 
U(l)+Higgs theory is of the first order Q. For 
larger x (i.e., for type II superconductors), the 
IR-problems are more severe since the transition 
gets weaker, and the order has remained unclear. 
Arguments in favour of a second order transition 
have been given, e.g., in but a still higher- 

order transition cannot at present be excluded. 
In particular, in the lattice studies in H||, the 
correlation lengths were not measured. 

It is interesting to compare the situation with 
that in the 3d SU(2)+Higgs theory. There the 
line of first-order transitions ends at x ~ 1/8, 
and for x > 1/8 the transition is of higher than 
second order jy]]. An important difference be- 
tween the U(l) and SU(2) cases is that in the lat- 
ter all the excitations in the symmetric phase are 
massive, whereas in U(l) there should always be 
a massless photon in the continuum limit ]6],[T^- 
|l~4| ]. Consequently, one can define an order pa- 
rameter for the continuum theory p"2[ . Away 
from the continuum limit an exponentially small 
non-perturbative mass is generated in the com- 
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Figure 2. Susceptibility at x = 0.0463. 



Figure 3. Susceptibility at x = 2. 



pact lattice formulatio n fl5|| , and the phases are 
analytically connected [|l6| . 

To study the order of the phase transition, we 
have used two different values of x. The first case 
x = 0.0463 corresponds to a strongly type I super- 
conductor, and the other case x = 2 to a strongly 
type II superconductor. 

In the analysis we follow closely jll]]. The 
method employed is to plot the maximum of the 
susceptibility as a function of the lattice volume. 
The susceptibility x 1S defined as 

x = e ^<(r^-(0*0)) 2 ). 

In a first order transition x grows as the volume 
V and in a second order transition the expected 
behaviour is ~ V K , n being a critical exponent. If 
k < 0, or if the transition is of higher than second 
order, x ~ V°. 

Our results are plotted in Fig. Jl, and the in- 
dividual susceptibilities in Figs. g-|3. The system 
has a first order phase transition at x = 0.0463, 
as the maximum of susceptibility grows linearly 
with volume. This is not the case with the x = 2 
data, so the transition in strongly type II super- 
condutors is not of first order. The data would 
suggest that if the transition is of second order 
the critical exponent k is close to zero, as noticed 
already in [pPJ]. However, the transition might 
also be of higher order. 



To study the transition at x = 2 in more de- 
tail, we have measured the masses of the scalar 
and vector excitations. The masses are measured 
from the exponential decay of suitably chosen cor- 
relators. For the scalar mass we used the operator 

<j>*(x)<j>(x), 

and for the vector mass 
Im <fi* (x)Ui(x)4>{x + i). 

The photon mass in the symmetric phase can 
be measured from the imaginary part of the pla- 
quette variable at non-zero momenta Jl7| . 

The scalar mass is shown in Fig. [| and the vec- 
tor mass in Fig. ^. The scalar mass is clearly 
non-vanishing. For the vector mass the signal 
gets very noisy in the symmetric phase {y > y c , 
@h < Ph,c), and the techniques of 18 should 
be used. In general, improved techniques make 
the masses smaller so that, for example, the vec- 
tor mass could go to zero at the phase transition 
point. Our very preliminary photon mass mea- 
surements in the symmetric phase are consistent 
with zero. 

In conclusion, we have seen that in type II su- 
perconductors the phase transition is not of first 
order. It may even not be of second order, since 
we see a non- vanishing scalar mass. For definite 
conclusions we will have to improve the quality of 
mass measurements and study the finite a effects. 
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Figure 4. Scalar mass near the pseudocritical cou- 
pling Ph,ci f° r x = 2. The values of /3h,c can be 
seen from the maxima of x m Fig- 3. 



4.0 



3.0 - 



.0^2.0 - 

E 



1.0 - 



0.0 



■32' 
1 48 s 




1. P. Arnold, Phys. Rev. D 46 (1992) 2628; K. 
Farakos, K. Kajantie, K. Rummukainen and 
M. Shaposhnikov, Nucl. Phys. B 425 (1994) 
67; A. Jakovac, A. Patkos and P. Petreczky, 
Phys. Lett. B 367 (1996) 283. 

2. B.I. Halperin, T.C. Lubensky and S.-K. Ma, 
Phys. Rev. Lett. 32 (1974) 292. 

3. C. Dasgupta and B.I. Halperin, Phys. Rev. 
Lett. 47 (1981) 1556. 

4. J. Bartholomew, Phys. Rev. B 28 (1983) 5378. 

5. J. March-Russell, Phys. Lett. B 296 (1992) 
364. 

6. W. Buchmiiller and O. Philipsen, Phys. Lett. 
B 354 (1995) 403. 

7. B. Bergerhoff, F. Freire, D.F. Litim, S. Lola 
and C. Wetterich, Phys. Rev. B 53 (1996) 
5734. 

8. K. Kajantie, M. Karjalainen, M. Laine and J. 
Peisa, in preparation. 

9. M. Karjalainen a nd J. Peisa, HU-TFT-96-27 
|hcp-lat/9607023| 1. 

10. M. Laine, Nucl. Phys. B 451 (1995) 484. 

11. K. Kajantie, M. Laine, K. Rummukainen and 
M. Shaposhnikov, CERN-TH/96-126 Jhep~ 



ph/960528S 



12. A. Kovner, B. Rosenstein and D. Eliezer, 
Nucl. Phys. B 350 (1991) 325. 

13. A. Hebecker, Z. Phys. C 60 (1993) 271. 

14. J. -P. Blaizot, E. Iancu and R.R. Parwani, 
Phys. Rev. D 52 (1995) 2543. 

15. A.M. Polyakov, Phys. Lett. B 59 (1975) 82. 

16. T. Banks and E. Rabinovici, Nucl. Phys. B 
160 (1979) 349; E. Fradkin and S. Shenker, 
Phys. Rev. D 19 (1979) 3682. 

17. B. Berg and C. Panagiotakopoulos, Phys. 
Rev. Lett. 52 (1984) 94; H.G. Evertz, K. 
Jansen, J. Jersak, C.B. Lang and T. Neuhaus, 
Nucl. Phys. B 285 (1987) 590. 

18. O. Philipsen, M. Teper and H. Wittig, Nucl. 
Phys. B 469 (1996) 445. 



Figure 5. Vector mass near /3h,c for x = 2. Note 
that the high-temperature symmetric phase cor- 
responds to /3 h < /3h,c- 



